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The adequate interpretation of scanning tunneling microscopy (STM) images of the clean Si(001) 
surface is presented. We have performed both STM observations and ab initio simulations of STM 
images for buckled dimers on the clean Si(001) surface. By comparing experimental results with 
theoretical ones, it is revealed that STM images depend on the sample bias and the tip-sample 
separation. This enables us to elucidate the relationship between the corrugation in STM images 
and the atomic structure of buckled dimers. Moreover, to elucidate these changes, we analyze details 
of the spatial distributions of the n, n* surface states and a, a* Si-Si bond states in the local density 
of states which contribute to STM images. 



I. INTRODUCTION 

Scanning tunneling microscopy (STM) has provided 
much useful information on surface configurations, with 
atomic resolution. The clean Si(001) surface is vigor- 
ously studied by STM because of the great fundamental 
and technological importance of semiconductor devices. 
In the STM images of the 2x1 structure on the terrace, 
silicon dimers appear to be symmetric at room temper- 
ature due to the time average of the flip-flop motion of 
buckled dimers, while at temperatures below ^200 K, 
the stable configuration is the c(4x2) oj-ai(2x2) structure 
consisting of asymmetric dimer rows.uu For the c(4x2) 
structure, the dependence of the STM image on the bias 
voltage at 80 K has been studied in detail.Q On the other 
hand, the behavior of the dimers at step edges differs 
from that on the terrace. For example, the dimers at Sa 
steps are observed to be buckled at both low and room 
temperatures!] However, detailed observation of the Sa 
steps has not been exhaustively carried out, since filled- 
and empty-state images at only a few bias voltages have 
been reported^ 

It is known that STM images reflect electronic struc- 
tures rather than atomic geometries. Consequently, bias- 
dependent STM images result from the variation of the 
spatial distribution of the local density of states (LDOS) 
which is a function of energy. Therefore, to interpret 
STM images properly, the spatial distributions of the 
LDOS must be calculated strictly by ab initio methods. 
Some authors have attempted to compare experimental 
STM images with theoretical ones for the Si(001) re- 
constructed surface. Kageshima and Tsukada carried 
out the simulation of polarity-dependent STM images 
of Si(001). However, their investigation was based on 
the semiFjempirical method with linearly combined atomic 
orbitals.Q Hata et al. pointed out that the tt, n* surface 
states and the a, a* bulk ones contribute to the STM 
images, by calculating the spatial distributions of those 
states localized near the outermost atomsfl These results 



emphasize that a reproduction of the STM image from 
the calculated spatial distribution of the LDOS is very 
important. But such research is still limited at present. 

In this paper, we present an improved interpretation of 
STM images and provide a more profound understand- 
ing of the electronic structure for the clean Si(001) sur- 
face. STM observations are performed on the buckled 
dimers at the Sa step of the clean Si(001) surface. They 
are strongly dependent on the applied bias voltage and 
the tip-sample separation, which is confirmed by ab ini- 
tio simulations. The STM images will be interpreted as 
isosurfaces of the spatial distribution of the LDOS in 
the vacuum region between a tip and a sample, which 
is greatly different from the behavior of the LDOS near 
the outermost atoms. Moreover, the calculated LDOS 
shows how the spatial distributions of the characteristic 
surface states (tt, ir^, and 7r| dangling-bond states) and 
bulk states (a and a* Si-Si bond states) contribute to the 
bias dependence of the STM images. 



II. METHODS OF STM OBSERVATION AND 
THEORETICAL CALCULATION 

A. STM observations 

Experiments were performed in an ultrahigh-vacuum 
(UHV) chamber equipped with STM. The base pressure 
was less than lxlCP 8 Pa. The samples were Si(001) 
(B-doped, 0.2-0.5 17cm) wafers. They were degassed at 
900 K overnight, flashed to 1450 K for 30 sec and then 
cooled to room temperature to prepare the clean Si(001)- 
2x1 surface. The tips were prepared by the electro- 
chemical etching of a polycrystalline tungsten wire and 
cleaned by electron bombardment in UHV prior to use. 
The constant-current STM images were obtained at room 
temperature. 
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B. Theoretical calculations 

The ab initio electronic-states calculations were per- 
formed by the self-consistent pseudoaotential method 
within the density-functional theory.13 The exchange- 
correlation interaction was treated by the Ceperley=Alder 
form in the local-density approximation (LDA)Ja The. 
norm-conserving Troullier-Martins pseudopotentialsL2l 
were used in a separable nonlocal formJlJ The wave func- 
tions were expanded in a plane- wave basis set with a cut- 
off energy of 24 Ry. 

The Si(001)-p(2x2) surface with asymmetric dimers 
was modeled as a repeated slab consisting of five sili- 
con layers, the lowest of which was terminated by hy- 
drogen atoms, and the width of the vacuum region was 
equivalent to that of 12 silicon layers. The four outer- 
most silicon layers were fully-optimized by the ab initio 
molecular-dynamics method.E3 

Many theoretical STM images have, been generated 
in the Tersoff-Hamann approximation.E3~t3 This method 
is actually—valid for many systems despite its extreme 
simplicity .Efl In this scheme, the tunneling current is pro- 
portional to the surface LDOS at the tip position inte- 
grated over the energy range restricted by the applied 
bias voltage. We calculated the LDOS p s (x,y, z;e) at 
spatial points (x, y, z) and energy e by sampling the 100 
fc-points of the surface Brillouin zone (SBZ). The STM 
images were generated from the isosurface of the spatial 
distribution formed by the integration of the LDOS over 
the energy range from Fermi energy Ep to Ep-eV with 
applied voltage V, i.e., J~ p s (x,y, z\ Ep — eV + e)de. 
Hereafter, we refer to this 'energy-integrated LDOS' as 
the 'EI-LDOS'. Thereby the computed STM images cor- 
respond to experimental constant-current images. 



III. RESULTS AND DISCUSSION 

A. Optimized atomic structure and electronic 
structure of the Si(001)-p(2x2) surface 

The atomic structure of the Si(001)-p(2x2) surface was 
optimized by minimizing the total energy. The dimer 
bond length of 2.35 A, which is the same length as the 
bond length in bulk silicon, and the buckling angle of 
19.0° were found. These values are similar to the dimer 
bond length of 2.38 A and the buckling angle of 18.4° for 
the p(2x2) reconstruction reported by Shkrebtii et al£3 
Our calculated dimer bondriength is longer than the 2.28 
Aiqund by Ramstad et alc3 and 2.31 A by Gunnella et 
alE3 The backbond lengths of the upper and lower atoms 
are 2.38 and 2.33 A, respectively, which are longen.than 
the of 2.34 and 2.31 A reported byJjlamstad et alt3 and 
2.33 and 2.30 A by Gunnella et aM 

We calculated the band structure for the optimized 
Si(001)-p(2x2) surface. This is formed at 100 equidistant 
points along the Bloch wave vector k|| in two-dimensional 



(001) SBZ. Figure [l] shows the band structure for the 
Si(001)-p(2x2) surface along high-symmetry directions 
of the SBZ. The single-particle calculation within the 
LDA is well known to underestimate the band gap. Al- 
though the band gap is calculated as 0.15 eV, all empty 
states of the band structure must shift rigidly upwards 
in energy by 0.65 eV.Eai23 This value was determined 
from the results of angle-resolved. photoemission£il angle- 
resolved inverse photejemissionJ23 and scanning tunneling 
spectroscopy (STS).cj The position of the Fermi energy 
(Ep) is determined at about 0.6 eV above the valence- 
band top and about 0.2 eV below the conduction-band 
bottom in order to coincide with our STM results. 

Filled 7r and empty tt* surface states separate into two 
filled (tti, 772) and empty (7r 1; tt^ ) surface states in the 
bulk band gap, respectively. This fea ture ., is consistent 
with the results of other calculations .ESE-aH The filled 
7Ti and 7T2 states disperse from -1.2 to -0.6 eV. Although 
the empty irl state widely spreads from +0.2 to +1.1 eV, 
the empty state spreads at around +1.1 eV. The 7Ti 
and 7T2 states, as well as the tt^ and tt^ states, overlap 
each other from the K to the J' point (in the direction 
perpendicular to the dimer rows). 



B. Comparison of experimental STM images with 
simulated ones 



We aimed at obtaining the STM images of a buckled- 
dimer row at the Sa step on the clean Si(001)-2xl sur- 
face under different sample-bias conditions. Figures ^ 
(filled states) and || (empty states) show the experimen- 
tal STM images of the Sa step and the simulated ones 
of the Si(001)-p(2 x 2). All STM images in Figs. | and | 
were obtained at a tunneling current of 0.5 nA, and sim- 
ulated images were generated by EI-LDOS of 8.2xl0 -5 
electrons/ A 3 . The simulated images accurately repro- 
duce the experimental results. 

In the experimental images of filled states [Figs. ||(a) 
and H(c)] , the dimer row at the Sa step makes a zigzag 
pattern. In other words, one atom in a buckled dimer is 
brightly observed and the other is absent £n3 The filled- 
state images do not change greatly in the range of bias 
voltage from -1.0 to -2.5 V, except for the difference in 
the corrugation height. The simulated filled-state images 
[Figs. |2|(b) and §(d)] imply that the bright protrusion 
of the buckled dimer in the experimental images corre- 
sponds to the upper atom of the buckled dimer. 

On the other hand, the empty-state images of the buck- 
led dimer at the Sa step show dramatic changes at room 
temperature as the sample bias is increased. These im- 
ages can be divided into three typical types depending 
on applied voltage: (i) below +1.1 V, (ii) at about +1.2 
V, and (iii) from +1.3 to +2.0 V. The empty-state im- 
ages in Fig. ^ were obtained for the same area, as shown 
in Fig. |H At a sample bias of +0.6 V in Fig. ||(a), the 
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zigzag corrugation reverses to, that in the filled-state im- 
ages in Figs. ||(a) and §(c).Q This reversal is also ob- 
served on the terrace for the Si(001)-c(4x2) surface at 
low temperaturesou Inside one buckled dimcr in Fig. [|, 
two protrusions are clearly distinguished. One is brighter 
than the other, which cannot be seen in the filled states 
[Figs. |(a) and |(c)]. The simulated STM image at a 
sample bias of +0.6 eV in Fig. ||(b) coincides with the 
experimental one [Fig. ^(a)]. The brighter protrusion in 
Fig. ||(a) corresponds to the lower atom of the buckled 
dimer. When +1.2 V is applied to the sample, the exper- 
imental and simulated images show that the zigzag dimer 
row changes into a symmetric one [Figs. 0(c) and |^(d)]. 
At a sample bias of +1.5 V in Fig. \^(e), the brighter 
oval-shaped protrusion is located at the opposite side to 
that at a sample bias of +0.6 V in the buckled dimer. 
The atomic protrusions align parallel to the dimer rows. 
In addition, a deep trough appears in the middle of the 
dimer row, as indicated by black arrows in Fig. [|(e) . This 
has been observed in the empty-state images, at high bi- 
ases at both room and low temperatures .BCj These fea- 
tures are accurately represented in the simulated image, 
as shown in Fig. ||(f). The bright protrusion corresponds 
to the upper atom of the buckled dimer. 

We also investigate the current dependence of buckled 
dimers and dimer rows on the Si(001) surface in STM 
images. Figure |(a) shows STM images of the Si(001)- 
2x1 surface at a sample bias of 1.5 V. In the middle of 
the image, the tunneling current was switched from 0.5 
to 5.0 nA. In the upper (lower) area of the image, the 
tunneling current was 0.5 (5.0) nA. In the upper half of 
the figure, the deep trough appears in the center of the 
dimer row, as mentioned above. The deep trough exists 
in the center of the dimer row indicated by black arrows 
in Fig. |](a) (at a tunneling current of 0.5 nA), as men- 
tioned above. When the tunneling current is changed 
to 5.0 nA, the other trough appears between the dimer 
rows, as indicated by white arrows in Fig. ||(a). Further- 
more, the two atoms in a dimer are clearly visible. Since 
the tip approaches the sample surface on switching the 
tunneling current from 0.5 to 5.0 nA, the LDOS at the 
nearer sample surface contributes to the tunneling cur- 
rent. In Fig. [|(c) the simulated image is generated as 
the isosurface of 8.2xl0 -4 electrons/A 3 , which is one or- 
der larger than that adopted in Fig. ||(f). This simulated 
image is consistent with the experimental one and indi- 
cates that STM images depend on the magnitude of the 
surface LDOS, i.e., the tunneling conductance. 

Consequently, the simulation of STM images can sup- 
port the interpretation of the experimental STM images, 
as shown in Figs. ||, ||, and || The simulated ones gen- 
erated by the EI-LDOS of 8.2xl0~ 5 electrons/A 3 cor- 
responded well with the experimental STM images ob- 
tained at the tunneling current of 0.5 nA. We can esti- 
mate the distance from the surface at which the spatial 
distribution of the EI-LDOS contributes to the tunneling 
current, i.e., the STM image. When the sample bias and 
the tunneling current are set to +1.5 and 0.5 (5.0) nA, 



the distance of the EI-LDOS from the surface contribut- 
ing to the STM images is predicted to be 4.6 (3.7) A. 
Similarly, when the tunneling current is fixed to be 0.5 
nA, the EI-LDOS at about 4.9, 4.7, 3.4, and 4.2 A from 
the surface contributes to the STM images at the sam- 
ple biases of -2.0, -1.0, +0.6, and +1.2 V, respectively. 
As the absolute value of the bias voltage is increased or 
the tunneling current is decreased, the distance from the 
surface at which the surface LDOS contributes to the 
tunneling current becomes large. When tip-sample sep- 
arations are assumed to be about twice as long as the 
distances estimated above, they agree with the experi- 
mental valuesE-X To estimate the tip-sample separation 
exactly, both the tip and the sample must be taken into 
account in the calculation. 



C. LDOS and spatial distribution of each state in 
LDOS of the Si(001)-p(2x2) surface 

By calculating the LDOS at different sites on the 
Si(001)-p(2x2) surface, we demonstrate how the spatial 
distributions of the characteristic surface and bulk states 
in the LDOS contribute to the STM images. 

Figure || shows the LDOS at different sites on the 
Si(001)-p(2x2) surface. The LDOS in Fig. |(b) are gen- 
erated at seven specific sites (sites A-G) in Fig. ||(a) 
which are located at the height of 4.3 A above the upper 
dimer atom. In the filled states, we find a sharp peak 
at -0.7 eV and a broad one at -1.7 eV marked by FSi 
and FS2 respectively, in Fig. |5|(b). These peaks appear 
similarly in every spectrum in Fig. ^(b). By comparing 
these peaks with the band structure (Fig.[l]), FSi and FS2 
mainly correspond to occupied [tt (tti and ^2)] dangling- 
bond states and the Si-Si bond (a) state, respectively. 
FSi and FS2 at the upper atom [site A in Fig. ||(a)] are 
the strongest, while they are the weakest at the cave and 
valley bridge sites [sites D and G in Fig. 0(a)] . 

For the empty states, the LDOS have a small peak at 
around +0.8 eV (ESi), a sharp one at +1.2 eV (ES 2 ), 
and a broad one at +1.8 eV (ES3), as shown in Fig. ||(b). 
From the band structure, ESi and ES2 mainly correspond 
to unoccupied irl and tt^ dangling-bond states, respec- 
tively. The heights and the fine structures of ESi and 
ES2 at each site greatly differ from one another, which 
causes a dramatic change in the empty-state STM im- 
ages. ESi at the lower atom, bridge, and pedestal sites 
[sites B, C, and E in Fig. ||(a)] are stronger than at other 
sites. ES2 at the upper atom, cave, side bridge, and valley 
bridge sites (sites A, D, F, and G) is much sharper than 
at other sites. In contrast, ES2 at the pedestal site (site 
E) decreases markedly. As described below, the analy- 
sis of the spatial distribution of each state is useful for 
understanding STM images. 

We compare our calculated LDOS with previ- 
ously reported— .STS measurements for the Si(001)- 
2x1 surfaceSEaO The experimental tunneling spectra 
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(dI/dV)/(I /V) have some characteristic features, (i) 
The band gap is ~0.6 eV. (ii) In the filled states, a large 
peak is found around -0.7 eV. (iii) In the empty states, 
two peaks are observed at around +0.6 and +1.3 eV. 
These features agree with those of our theoretical LDOS. 
We attribute the peaks at -0.7, +0.6, and +1.3 eV to n, 
7r* , and , respectively. 

We calculate the spatial distribution of each character- 
istic state by integrating the LDOS over the energy range 
in which peaks relating to the state spread are found. We 
also define this 'energy-integrated LDOS' as 'EI-LDOS'. 
Figure || shows the spatial distributions of EI-LDOS of 
the two characteristic peaks [(a) FSi and (b) FS2] in the 
filled states. Figure. ||(a) shows that the charge transfer 
occurs from the lower atom to the upper one aad the 
EI-LDOS is largely localized at the upper atomE3. The 
maximum EI-LDOS in the vacuum region is located just 
above the upper atom. For FS2 [Fig. ||(b)], the spatial 
distribution of the EI-LDOS is largely localized between 
Si-Si atoms and forms a bonds, which is greatly differ- 
ent from the tt surface state (FSi). In contrast, in the 
vacuum region, the maximum EI-LDOS of FS2 is again 
located above the upper atom, although the corrugation 
of contour lines for FS2 is smaller than that for FSi . This 
feature leads to the similarity of STM images at any sam- 
ple bias in the filled state, as shown in Fig. |^. 

Figure [?] shows the spatial distributions of EI-LDOS of 
the three peaks [(a) ESi, (b) ES2, and (c) ES3] in the 
empty states. In Fig. g(a), the EI-LDOS of ESi is lo- 
calized at the lower atom of the buckled dimer. In the 
vacuum region, the corrugation of the EI-LDOS has two 
humps above a buckled dimer and a minimum at the 
center of the dimer. The hump is considerably larger 
above the lower atom than above the upper one. The 
ir"l contributes to the STM images at low biases (below 
+1.1 eV). On the other hand, the spatial distribution of 
ES2 is extremely different from that of ESi. As shown 
in Fig. ^(b), the deep minimum appears in the middle 
of a buckled dimer [sites C and E in Fig. ||(a)] and the 
maximum of the EI-LDOS of ES2 extends into the vac- 
uum above the upper atom more than above the lower 
atom. These striking features of the tt^ state lead to the 
deep trough in the middle of a dimer row and the bright 
oval-shaped protrusion at the upper dimer atom in the 
STM images at the sample bias from +1.2 to +2.0 V, as 
shown in Figs. ||(e) and 0(f). In addition, we also find 
that the EI-LDOS of ES2 decays slower in vacuum at 
the cave, valley bridge, and side bridge sites (sites D, F, 
and G) than at the other sites. This feature explains the 
changes in the STM image in Fig. ||(a) when a tunneling 
current is switched from 0.5 to 5.0 nA at a sample bias 
of +1.5 V. In other words, since the tip-sample separa- 
tion is shorter and the STM image reflects the spatial 
distribution of the LDOS nearer to the surface, the other 
trough emerges between dimer rows and each atom ap- 
pears to be clearly resolved. ES3 corresponds to the a* 
antibonding bulk state. As shown in Fig. f7|(c), in the 
bulk, this state distributes not only at dimer rows but 



also between dimer rows and generally does not local- 
ize near the dimer atoms, in contrast with the n\ and 
7i"2 states. In the vacuum region, the corrugation of the 
EI-LDOS extends slightly above the upper atom. When 
the sample bias is higher, the position of the spatial dis- 
tribution forming the STM image becomes farther from 
the surface. The LDOS of the a* state decays slower in 
vacuum than that of the ~k\ and it^ states. Therefore, we 
predict that the a* state is evident in the STM images 
at biases higher than +2.1 V. 

We have presented the features of the spatial distribu- 
tion of each state in the LDOS and revealed how the sur- 
face and bulk states contribute to STM images. We also 
demonstrated that the STM image of the clean surface 
can be interpreted as the isosurface of the spatial distri- 
bution of the EI-LDOS in the vacuum region, which is 
greatly different from that near the outermost atoms. 



IV. CONCLUSION 

We have performed STM observations and theoretical 
simulations of STM images for the buckled dimer on the 
clean Si(001) surface. By comparing the experimental re- 
sults with theoretical ones, the relationship between STM 
images and the atomic structure of buckled dimers of the 
Si(001) surface was clarified. The spatial distributions of 
the surface and Si-Si bond states in the LDOS were in- 
vestigated. For the filled states, the ir (m and 7^) surface 
states contribute to the low-bias STM images and, as the 
sample bias is increased, the Si-Si bond (ex) state is added 
to the 7r states. Since the corrugations of the spatial dis- 
tributions of the 7r and a states are similar in vacuum, the 
filled-state images are nearly identical to each other. On 
the other hand, the empty-state images change consider- 
ably with relative to the sample bias. The STM images 
obtained at low biases reflect the spatial distribution of 
the 7T* surface states. As the sample bias is increased, the 
7T2 surface state begins to dominate the STM images, and 
a deep trough appears at the center of the dimer at high 
biases above +1.3 eV. Moreover, when the tunneling cur- 
rent is set to be larger, the LDOS nearer to the sample 
surface contributes to the STM image. These analyses 
enable us to understand not only STM images but also 
the electronic structure of the clean Si(001) surface. 
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K J 
Wave Vector 

FIG. 1. Band structure of the Si(001)-p(2x2) surface op- 
timized by molecular dynamics. Eigenvalues in the empty 
states are rigidly shifted by 0.65 eV. To match our STM re- 
sults, the Fermi energy Ef position is determined at 0.6 eV 
above the valence band top. 
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Vs = -2.0 V 




Vs = -1 .0 V 




FIG. 2. Filled-state STM images (5x5 nm 2 ) of the Sa step 
on the Si(001)-2x 1 surface at room temperature (left column) 
and simulated ones of the Si(001)-p(2x2) surface (right col- 
umn). These experimental images were obtained at sample 
biases of (a) -2.0 and (c) -1.0 V and at a tunneling current 
of 0.5 nA. These simulated images were generated from the 
isosurfaces of EI-LDOS [(b) to -2.0 and (d) to -1.0 eV] 
of 8.2xl0 -5 electrons/A 3 . The maximal heights of these iso- 
surfaces are (b) 4.9 and (d) 4.7 A from the upper atom. The 
atomic positions in these images are represented in (e) and 
the smallest filled circles are the second-layer atoms. 




Vs = +0.6 V 




Vs = +1 .2 V 




Vs = +1 .5 V 

FIG. 3. Empty-state STM images (5x5 nm 2 ) of the Sa 
step on the Si(001)-2xl surface at room temperature (left col- 
umn) and simulated ones of the Si(001)-p(2x2) surface (right 
column). These experimental images are obtained at sample 
biases of (a) +0.6, (c) +1.2, and (e) +1.5 V and at a tunneling 
current of 0.5 nA. These simulated images are generated at 
the isosurfaces of EI-LDOS [(b) to +0.6, (d) to +1.2, and 
(f) to +1.5 eV] of 8.2xl0~ 5 electrons/ A 3 . The maximal 
heights of these isosurfaces are (b) 3.4, (d) 4.2, and (f) 4.6 A 
from the upper atom. The atomic positions in these images 
are represented in Fig. ^|(e). 




8.2x 10-5 ( e /A3) 8.2x 10-4( e /A3) 

FIG. 4. Current-dependent STM images (9x9 nm 2 ) of the 
Si(001)-2xl surface of the empty states at room temperature 
and the simulated image, (a) The experimental image is ob- 
tained at a sample bias of +1.5 V. In the upper and lower parts 
of the STM image, the tunneling current is 0.5 and 5.0 nA, re- 
spectively, (b) This image is equivalent to Fig. |^(f). (c) The 
simulated image is generated by the isosurface of 8.2 xlO -4 
electrons/ A 3 , which is one order larger than that adopted in 
Fig. |^(f). The maximal height of this isosurface is 3.7 A from 
the upper atom. 




-2-1012 
Energy (eV) 

FIG. 5. Calculated LDOS of the Si(001)-p(2x2) surface at 
different sites, (a) Top view indicates the seven sites at which 
the LDOS are calculated. A, upper atom; B, lower atom; C, 
bridge; D, cave; E, pedestal; F, side bridge; and G, valley 
bridge, (b) The LDOS at these seven sites are generated at 
the same height of 4.3 A from the upper atom. The five char- 
acteristic peaks at -0.7, -1.7, +0.8, +1.2, and +1.8 eV are 
marked by FSi, FS2, ESi, ES2, and ES3, respectively. 
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FIG. 6. Isosurfaces (left panel) and contour plots (right 
panel) of charge densities for the two characteristic peaks in 
the filled states: (a) -0.7 (-0.6 to -1.3) and (b) -1.7 (-1.4 to 
-2.0) eV. Isosurfaces (left panel) in the bulk (white) and the 
vacuum (at ~4.3 A from the upper atom; green) are shown in 
the same cell; the height of silicon layers is represented by the 
size of the ball. Logarithmic contour maps (right panel) are 
along the cross section in the (110) plane along the [110] direc- 
tion including the buckled dimer. The magnitude of the en- 
ergy-integrated LDOS increases/decreases one order per four 
contour lines. 




(») ES< x ,' ft)* {< 2 19 +1 *V 




FIG. 7. Isosurfaces (left panel) and contour plots (right 
panel) of charge densities for the three characteristic peaks 
in the empty states: (a) +0.8 (+0.2-+1.0), (b) +1.2 
(+1.1— 1-1.3), and (c) +1.8 (+1.4-+2.0) eV. The meanings 
of the symbols are the same as in Fig. H. 
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